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MILD  (flameless)  combustion  is  a  new  rapidly  developing  technology.  The  IFRF  trials  have  demonstrated 
high  potential  of  this  technology  also  for  N-containing  fuels.  In  this  work  the  IFRF  experiments  are  ana¬ 
lyzed  using  the  CFD-based  mathematical  model.  Both  the  Chemical  Percolation  Devolatilization  (CPD) 
model  and  the  char  combustion  intrinsic  reactivity  model  have  been  adapted  to  Guasare  coal  combusted. 
The  flow-field  as  well  as  the  temperature  and  the  oxygen  fields  have  been  accurately  predicted  by  the 
CFD-based  model.  The  predicted  temperature  and  gas  composition  fields  have  been  uniform  demonstrat¬ 
ing  that  slow  combustion  occurs  in  the  entire  furnace  volume.  The  CFD-based  predictions  have  high¬ 
lighted  the  NOx  reduction  potential  of  MILD  combustion  through  the  following  mechanism.  Before  the 
coal  devolatilization  proceeds,  the  coal  jet  entrains  a  substantial  amount  of  flue  gas  so  that  its  oxygen 
content  is  typically  not  higher  than  3-5%.  The  volatiles  are  given  off  in  a  highly  sub-stoichiometric  envi¬ 
ronment  and  their  N-containing  species  are  preferentially  converted  to  molecular  nitrogen  rather  than  to 
NO.  Furthermore,  there  exists  a  strong  NO-reburning  mechanism  within  the  fuel  jet  and  in  the  air  jet 
downstream  of  the  position  where  these  two  jets  merge.  In  other  words,  less  NO  is  formed  from  combus¬ 
tion  of  volatiles  and  stronger  NO-reburning  mechanisms  exist  in  the  MILD  combustion  if  compared  to 
conventional  coal  combustion  technology. 

©  2009  Published  by  Elsevier  Inc.  on  behalf  of  The  Combustion  Institute. 


1.  Introduction  and  objectives 

Coal  is  an  abundant  fuel  resource  in  many  of  the  developing 
regions  and  forecasts  show  that  it  is  likely  to  remain  a  dominant 
fuel  for  electricity  generation  in  many  countries  for  years  to  come. 
Coal-fired  power  plants  currently  generate  approximately  40%  of 
the  world  electricity.  Since  coal  dominates  the  energy  supply  in 
the  developing  countries  and  still  is  an  important  fuel  in  the  indus¬ 
trialized  nations  it  will  continue  to  play  an  important  role  in 
worldwide  power  generation  [1].  Thus,  the  development  of 
advanced  coal  combustion  technologies  of  a  higher  performance 
efficiency  and  lower  pollutant  emissions  is  a  major  goal  of  combus¬ 
tion  researchers. 

MILD  (Moderate  and  Intensive  Low-oxygen  Dilution)  combus¬ 
tion  (often  written  simple  as  mild  combustion),  named  also 
FLameless  Oxidation  -  FLOX  or  High  Temperature  Air  Combustion  - 
HTAC,  is  probably  the  most  important  achievement  of  combustion 
technology  in  recent  years.  Within  a  decade  or  two,  it  has  been 
developed  from  laboratory  tests  to  industrial  applications  which 
certainly  is  an  extraordinary  progress  as  for  an  energy  technology. 
A  comprehensive  review  paper  summarizing  both  the  develop¬ 
ment  and  current  status  of  this  technology  does  not  exist  although 
some  information  can  be  found  in  [2-5].  The  essence  of  this  tech¬ 
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nology  is  that  fuel  is  oxidized  in  an  environment  that  contains  a 
substantial  amount  of  inert  (flue)  gases.  For  clean  gaseous  fuels 
that  do  not  contain  any  fuel-bounded  nitrogen  this  results  in  very 
low  NOx  emissions  even  if  the  combustion  air  stream  is  preheated 
to  temperatures  in  excess  of  1000  °C.  Chemical  reactions  take  place 
in  almost  the  entire  volume  of  the  combustion  chamber  and  uni¬ 
formity  of  both  the  temperature  and  the  chemical  species  fields 
are  characteristics  of  this  technology. 

2.2.  MILD  combustion  of  coal 

The  International  Flame  Research  Foundation  (IFRF)  carried  out 
perhaps  the  first  experiment  on  applicability  of  the  MILD  combus¬ 
tion  concept  to  pulverized  coal  firing  [6-8].  A  high  volatile  coal 
(1.4%  N)  was  combusted  with  high  temperature  air  and  the  lowest 
NOx  emissions  were  in  the  range  160-175  ppm  (at  3%  02)  indicat¬ 
ing  high  NOx  reduction  potential  of  this  technology  also  for  nitro¬ 
gen  containing  fuels  (for  details  see  paragraph  2  below). 

A  pulverized  coal  FLOX  burner  is  currently  under  development 
in  the  frame  of  an  EU  project  [9].  The  tests  have  been  carried  out  in 
a  drop  tube  furnace,  in  a  bench-scale  reactor  (5-8  kW)  at  the  Uni¬ 
versity  of  Stuttgart  and  in  a  pilot-scale  rig  (40-60  kW)  at  RWTH 
Aachen  [10,11].  The  burner  combines  the  FLOX  technique  together 
with  air  staging.  The  primary  air  and  pulverized  coal  are  injected 
into  the  furnace  via  a  central  pipe.  The  secondary  air  is  provided 
through  nozzles  located  on  the  circumference  of  the  primary  jet 
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while  the  tertiary  (burnout)  air  is  injected  downstream.  Under 
optimized  firing  conditions,  the  NO*  concentrations  were  below 
200  mg/m3  (at  6%  02)  and  the  carbon  monoxide  concentrations 
were  as  low  as  30  mg/m3  (at  6%  02).  While  firing  Lausitz  lignite  a 
high  degree  of  burnout  (carbon  in  ash  <0.3%)  has  been  measured, 
however  burnout  problems  have  occurred  for  hard  Budryk  and  Kle- 
inkopje  coals  (carbon  in  ash  <18.9%)  [10,11].  Coal  behavior  under 
MILD  combustion  conditions  has  also  being  researched  at  the 
Ruhr-University  Bochum  [12,13]  to  gain  a  better  understanding 
of  NOx  formation  and  reduction  mechanisms.  The  data  on  volatile 
matter  release,  volatiles  composition,  and  char  are  determined  by 
flash  pyrolysis  investigations. 

MILD  combustion  of  pulverized  coal  was  experimentally  inves¬ 
tigated  by  Suda  et  al.  [14].  Flame  stability,  ignition  delay,  char 
burnout,  and  NOx  emission  were  examined  using  a  250  kW  fur¬ 
nace.  The  combustion  air  was  preheated  to  a  temperature  of  either 
623  or  1073  K.  The  authors  concluded  that  the  ignition  delay  de¬ 
creased  dramatically  with  increasing  air  temperature.  Stable  com¬ 
bustion  was  obtained  even  for  an  anthracite.  The  peak  temperature 
of  the  flames  remained  unaltered  when  the  air  preheat  was  in¬ 
creased  from  623  to  1073  K.  NOx  emissions  of  1000  and  400  ppm 
were  measured  for  the  623  and  1073  K  air  temperature,  respec¬ 
tively.  Thus,  a  40%  NOx  emissions  reduction  was  observed  when 
the  combustion  air  temperature  increased.  However,  the  absolute 
NOx  emissions  were  substantially  larger  than  these  obtained  in 
the  IFRF  experiment  [6]. 

Numerical  studies  of  MILD  coal  combustion  were  performed  by 
He  et  al.  [15]  who  modeled  Suda’s  experiments  [14].  The  goal  was 
to  study  NOx  formation  and  destruction  mechanisms.  The  free  jet 
theory  was  used  to  describe  the  flow  field.  It  was  assumed  that 
during  pyrolysis  the  coal  nitrogen  was  released  to  volatiles  as 
HCN  or  remained  in  the  char.  The  thermal  mechanism  of  NOx  for¬ 
mation  was  also  considered.  The  model  incorporated  three  NO 
destruction  paths:  the  reactions  with  HCN,  char  and  soot.  It  was 
concluded  that  90%  of  the  NOx  emissions  were  formed  through 
the  fuel  path  and  the  remaining  10%  through  the  thermal  route. 

The  experiment  on  MILD  coal  combustion  under  high  pressure 
(3  bar)  conditions  was  carried  out  by  Heil  et  al.  [16]  and  modeled 
by  Erfurth  et  al.  [17]  using  the  CFD  FLUENT  code.  Three-dimen¬ 
sional  steady-state  simulations  of  a  1  / 6-sector  of  the  furnace  were 
performed  for  lignites  and  bituminous  coals.  Standard  sub-models 
implemented  in  the  FLUENT  code  were  used;  the  Eddy  Dissipation 
Concept  with  two  global  reactions  for  turbulence-chemistry  inter¬ 
action  modeling  and  PI  or  Discrete  Ordinates  models  for  radiation. 
A  Lagrangian  description  for  the  solid  phase  was  used.  Very  simple 
empirical  sub-models  were  used  for  devolatilization  (constant 
rates)  and  char  burnout  (diffusion-kinetics  limited).  This  simple 
mathematical  model  was  able  to  predict  well  the  flow  field  and 
the  recirculation  inside  the  combustion  chamber.  The  tempera¬ 
tures  were  over-predicted  in  comparison  with  the  experimental 
data  while  the  species  concentrations  differed  substantially  from 
the  measured  values.  Needs  for  detailed  sub-models  for  devolatil¬ 
ization  and  char  burnout  became  apparent. 

The  IFRF  experiments  on  MILD  coal  combustion  technology  [6] 
were  simulated  by  Kim  et  al.  [18,19].  The  objective  was  to  investi¬ 
gate  the  ability  of  the  CFD  AIOLOS  code  to  predict  the  nitrogen  oxide 
levels.  Three-dimensional  simulations  of  a  quarter  of  the  furnace 
were  performed  under  steady-state  conditions.  The  coal  combus¬ 
tion  was  modeled  including  primary  and  secondary  pyrolysis,  vola¬ 
tile  combustion,  and  char  burnout,  using  sub-models  described  by 
Fortsch  [20,21].  All  the  gas  phase  reactions  were  calculated  using 
global  mechanisms.  The  pyrolysis  gas  was  represented  by  CH4  while 
the  tar  was  approximated  as  CxHyOz.  The  secondary  reactions  of  tar 
were  modeled  by  three  competing  reactions  yielding  soot,  hydro¬ 
gen,  light  hydrocarbons,  and  carbon  monoxide.  The  pyrolysis  prod¬ 
ucts  were  predicted  using  the  Functional  Group-Depolymerization 


Vaporization  Cross  linking  (FG-DVC)  model  [22]  which  allowed  for 
predictions  of  both  the  yield  and  the  composition  of  products  of  pri¬ 
mary  pyrolysis.  The  NOx  formation  modeling  included  fuel  and  ther¬ 
mal  mechanisms.  The  release  of  volatile  nitrogen  was  assumed  to 
follow  the  rate  of  primary  pyrolysis  and  HCN  and  NH,  were  the  pri¬ 
mary  N-containing  species  of  the  pyrolysis  gas.  The  release  of  char 
nitrogen  was  enhanced  in  the  initial  stages  of  char  combustion 
and  char  nitrogen  was  released  not  only  as  intermediate  NH,  but 
also  directly  as  NO,  depending  on  a  mechanism  factor  defined  by 
Fortsch  [20].  The  char-NO  production  was  proportional  to  the  rate 
of  CO  formation  in  char  oxidation.  The  interaction  between  chemis¬ 
try  and  turbulence  was  modeled  using  the  Eddy  Dissipation  Con¬ 
cept.  Turbulence  was  calculated  using  the  standard  k-e  model 
whilst  radiative  heat  transfer  by  a  Discrete  Ordinates  method.  An 
Eulerian  description  for  the  solid  phase  was  used.  Generally,  good 
agreement  between  experimental  data  and  simulation  results  was 
achieved.  However,  the  carbon  monoxide  concentrations  were 
strongly  over-estimated  inside  the  combustion  chamber  and  at 
the  furnace  exit.  Throughout  the  furnace  and  at  the  furnace  exit 
too  high  oxygen  concentrations  were  computed  while  the  carbon 
dioxide  concentrations  were  under-estimated.  The  numerical  mod¬ 
el  slightly  under-estimated  the  NOx  values  within  the  combustion 
zone.  The  authors  concluded  that  there  was  a  need  for  the  develop¬ 
ment  of  detailed  coal  combustion  models. 

2.2.  Objectives 

The  experiments  [6-8]  resulted  not  only  in  the  lowest  NOx 
emissions  achieved  in  the  40-years  long  IFRF  history  of  burner 
development  for  pulverized  coal  combustion  but  also  initiated  a 
new  burner/furnace  design  methodology.  Thus,  it  is  important  to 
gain  a  good  insight  into  the  NOx  formation  and  destruction  mech¬ 
anisms  of  the  MILD  combustion  in  order  to  identify  features  that 
distinct  this  new  technology  from  conventional  low-NOx  firing. 
To  obtain  such  an  insight  we  use  our  CFD-based  mathematical 
modeling  software  and  the  IFRF  experiments  [6-8,23].  At  the  same 
time,  a  comparison  of  the  model  predictions  against  the  measured 
data  provides  a  validation  of  several  sub-models  as  for  their  pre¬ 
dictability  of  this  emerging  technology. 

2.  The  IFRF  experiment 

The  experiments  [6-8]  were  executed  in  a  refractory  lined  IFRF 
furnace  No.  1  shown  in  Fig.  1.  The  burner  consisted  of  a  central 
125  mm  diameter  pipe  supplying  the  oxidizer  and  two  27.3  mm 
coal  injectors  (pipes)  which  were  located  280  mm  away  from  the 
burner  center,  as  shown  in  Fig.  1.  The  burner  operated  at 
0.58  MW  fuel  input  and  the  vitiated  air  temperature  was  1623  K, 
as  shown  in  Table  1.  The  vitiated  air  stream  contained  22%  (wet, 
by  weight)  oxygen  and  89  ppm  (wet,  by  weight)  NO.  The  vitiated 
air  and  the  coal  streams  were  supplied  into  the  furnace  with  injec¬ 
tion  velocities  of  65  and  26  m/s,  respectively. 

The  measurements  were  taken  at  several  traverses  (Fig.  1)  in  a 
horizontal  plane  cutting  through  the  burner  centerline.  In-flame 
measurements  included  mean  and  rms  axial  velocities,  gas  tem¬ 
perature,  gas  composition  (C02,  02,  CO,  NOx,  CnHm),  burnout,  solids 
concentrations,  total  radiance  and  total  radiative  fluxes  at  the  fur¬ 
nace  wall.  The  velocity  measurements  were  performed  using  the 
IFRF  water  cooled  Laser  Doppler  Velocimetry  (LDV)  probe.  In-flame 
temperatures  were  measured  using  a  suction  pyrometer  equipped 
with  a  type  B  thermocouple  (Pt6%Rh/Pt30%Rh).  Local  in-flame  gas 
compositions  were  measured  using  a  gas  sampling  probe.  The 
measurements  of  total  radiance  and  total  radiative  heat  fluxes 
were  performed  using  a  standard  IFRF  narrow  angle  radiometer 
and  an  ellipsoidal  radiometer,  respectively.  Detailed  description 
of  these  measurements  can  be  found  in  [6,7]. 
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Fig.  1.  Geometry  of  the  IFRF  furnace  No.  1;  location  of  the  measurement  traverses. 


Table  1 

Experimental  conditions. 


Mass  flow  (kg/h) 

Vel.  (m/s) 

Temp.  (K) 

Enthalpy  (MW) 

Composition  (wt%  wet) 

Coal 

66 

— 

— 

0.58 

Coal  transport  air 

130 

26 

313 

- 

02  =  23,  N2  =  77 

Oxidizer  (vitiated  air) 

675 

65 

1623 

0.30 

02  =  22,  H20  =  9.5,  C02 

=  12.5,  NO  =  89  x  10  4,  N2  =  56 

Furnace  exit 

871 

30 

1503 

0.36 

02  =  2.8,  H20  =  11,C02 

=  31,  NO  =  250  x  10~4,  N2  =  55.2 

3.  Mathematical  model 

The  CFD  FLUENT  code  (version  6.3)  has  been  used  in  our  work 
as  the  modeling  software.  Substantial  efforts  have  been  allocated 
into  a  proper  numerical  description  of  the  coal  combustion.  To  this 
end  several  fuel  specific  subroutines  have  been  developed  and 
these  are  described  in  detail  in  Section  3.1.  Sections  3.2-3.4  pro¬ 
vide  just  a  very  brief  description  of  the  sub-models  used. 


3.1.  Numerical  description  of  the  coal  combustion 

The  Venezuelan  bituminous  Guasare  coal  was  combusted  [6]. 
According  to  the  ASTM  classifications  it  is  a  high  volatile  bitumi¬ 
nous  A  coal.  The  proximate  and  ultimate  analysis  of  Guasare  coal 
are  given  in  Tables  2  and  3,  respectively. 

The  coal  was  milled  to  give  a  particle  size  distribution  of 
80%  <  90  pm.  The  Rosin-Rammler  distribution  function  with  the 
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Table  2 

Guasare  coal  proximate  analysis. 


Composition 

wt% 

Moisture  (105  °C) 

2.9 

Volatile  matter 

37.1 

Fixed  carbon 

56.7 

Ash 

3.3 

LCV 

31.74  (MJ/kg) 

Table  3 

Guasare  coal  ultimate  analysis  (dry, 

ash  free  basis). 

Composition  (wt%  daf) 

Coal 

Char 

Volatiles 

C 

81.6 

92.6 

72.51 

H 

5.5 

1.3 

9.10 

N 

1.5 

1.7 

1.3 

O 

10.7 

4.0 

16.3 

S 

0.6 

0.4 

0.8 

mean  diameter  of  42  jam  and  the  spread  of  1.36  represents  nicely 
the  measured  data,  as  shown  in  Fig.  2.  In  the  CFD  predictions  of  this 
work,  the  10-300  pm  size  range  has  been  divided  into  20  size  clas¬ 
ses  and  around  750  particles  are  injected  per  class  each  time  the 
particle  tracking  procedure  is  activated. 

In  contrast  to  the  empirical  devolatilization  models  (for  exam¬ 
ple  a  single  rate  model  [24],  Kobayashi  model  [25]),  the  Chemical 
Percolation  Devolatilization  (CPD)  model  [26-28]  characterizes 
the  devolatilization  behavior  of  a  rapidly  heated  coal  particle  on 
the  basis  of  the  physical  and  chemical  transformations  of  the  struc¬ 
ture.  Modeling  the  cleavage  of  the  bridges  and  the  generation  of 
light  gas,  char,  and  tar  precursors  is  then  considered  to  be  analo¬ 
gous  to  the  chemical  reaction  scheme  shown  in  Fig.  3. 

The  variable  £  represents  the  original  population  of  labile 
bridges  in  the  coal  lattice.  Upon  heating,  these  bridges  become 
the  set  of  reactive  bridges.  For  the  reactive  bridges  £*,  two  compet¬ 
ing  paths  are  available.  In  one  path,  the  bridges  react  to  form  side 
chains,  3.  The  side  chains  may  detach  from  the  aromatic  clusters  to 
form  light  gas,  gA.  As  bridges  between  neighboring  aromatic  clus¬ 
ters  are  cleaved,  a  certain  fraction  of  the  coal  becomes  detached 
from  the  coal  lattice.  These  detached  aromatic  clusters  are  the  hea- 
vy-molecular-weight  tar  precursors  that  form  the  metaplast.  The 
metaplast  vaporizes  to  form  coal  tar.  While  waiting  for  vaporiza¬ 
tion,  the  metaplast  can  also  reattach  to  the  coal  lattice  matrix 
(crosslinking).  In  the  other  path,  the  bridges  react  and  become  a 
char  bridge,  c,  with  the  release  of  an  associated  light  gas  product, 
g2.  The  total  population  of  bridges  in  the  coal  lattice  matrix  can 
be  represented  by  the  variable  p,  where  p  =  £  +  c. 


Fig.  3.  Scheme  of  the  coal  behavior  during  devolatilization  process  as  a  simplified 
network  of  chemical  bridges  [27]. 


Table  4 

The  parameters  for  the  CPD  devolatilization  model  of  Guasare  coal. 


Parameter 

Symbol 

Value 

Unit 

Initial  fraction  of  bridges  in  coal  lattice 

Po 

0.5 

— 

Initial  fraction  of  char  bridges 

c0 

0 

- 

Lattice  coordination  number 

0  + 1 

5 

- 

Cluster  molecular  weight 

Md 

300 

kg/kmol 

Side  chain  molecular  weight 

Mdel 

30 

kg/kmol 

There  are  five  coal-specific  parameters  of  the  CPD  model:  the 
average  molecular  weight  per  aromatic  cluster  (Md),  the  average 
molecular  weight  per  side  chain  (Mdd),  the  average  number  of 
attachments  per  cluster,  referred  to  as  the  coordination  number 
(<7  +  1 ),  the  fraction  of  attachments  that  are  bridges  (p0)  and  initial 
fraction  of  char  bridges  (c0),  as  shown  in  Table  4.  The  question  is 
how  to  determine  these  five  parameters  for  the  considered  coal? 
To  this  end  we  use  other  IFRF  experimental  data  [23]  on  character¬ 
ization  of  both  devolatilization  and  char  combustion  of  Guasare 
coal.  Fig.  4(left)  shows  the  devolatilization  curve  obtained  using 
the  IFRF  Isothermal  Plug  Flow  Reactor  (IPFR)  operated  at  a 
1200  °C  temperature.  The  particle  heating  rate  of  these  experi¬ 
ments  is  in  the  range  105  -f  106  K/s  depending  on  the  particle  size. 
Under  such  rapid  heating  conditions  a  55-75%  high  temperature 
yield  is  measured  as  opposed  to  37%  ASTM  volatiles.  The  CPD  dev¬ 
olatilization  model  has  been  run  using  as  input  the  eight  experi¬ 
mentally  determined  yields,  shown  in  Fig.  4(left).  The  initial 
fraction  of  char  bridges  has  been  fixed  to  zero  which  is  a  typical  va¬ 
lue  for  bituminous  coals.  The  obtained  values  of  the  other  CPD 
parameters  are  listed  in  Table  4  while  the  CPD  model  devolatiliza¬ 
tion  curve  is  shown  in  Fig.  4(left)  using  a  solid  line.  The  volatiles 
are  represented  in  the  CFD  predictions  as  Q .20H4.4sO0.44  which  pro¬ 
vides  the  stoichiometric  coefficient  of  1.5  (see  Eq.  (12)). 


parameter 

value 

unit 

mean  diameter 

42 

tim 

max.  diameter 

300 

tim 

min.  diameter 

10 

tim 

spread  parameter 

1.36 

Fig.  2.  Guasare  coal  particle  distribution  (left)  and  distribution  parameters  (right). 
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- model 

+  measurements 


- model  1400°C 

- model  1300°C 

model  950°C 

+  measurements  1400°C 
+  measurements  1300°C 
measurements  950°C 


Fig.  4.  Devolatilization  and  burnout  measurements  [23]  with  the  CPD  (left)  and  the  intrinsic  (right)  model  fittings  for  Guasare  coal. 


Char  particles  produced  from  rapid  pyrolysis  are  microporous 
solids  whose  properties  can  be  described  by  their  size,  true  and 
apparent  density,  porosity,  pore  volume  distribution  and  surface 
area  distribution.  The  rate  of  char  oxidation  is  controlled  by 
sequential  or  parallel  processes  of  oxygen  boundary  layer  diffusion, 
chemical  reaction  and  pore  diffusion.  The  intrinsic  model  for  char 
burnout  used  in  this  work  is  based  on  the  Smith’s  [29]  macroscopic 
pore  model  which  assumed  that  the  char  oxidation  reaction: 


with  pp  and  pt  being,  respectively,  the  apparent  and  true  char  den¬ 
sities  and  t  is  the  tortuosity  of  the  pores.  A  typical  value  for  t  is  V2 
which  corresponds  to  an  average  intersecting  angle  between  the 
pores  and  the  external  surface  of  45°.  The  Knudsen  diffusion  coeffi¬ 
cient  D/cn  is  calculated  as 


Dkr  =  970 rp 


C(s)  +  02(g)  =  C02(g) 


is  of  the  first  order.  The  overall  surface  reaction  rate  includes  the  ef¬ 
fect  of  bulk  diffusion  and  chemical  reaction.  The  diffusion  coeffi¬ 
cient  D0  is  computed  via: 


Do  —  Ci 


(Tp  +  T^/2 

dP 


The  chemical  rate  is  expressed  in  terms  of  the  intrinsic  chemical 
and  pore  diffusion  rates: 


where  dp  is  the  particle  diameter,  r\  the  effectiveness  factor,  Tp  the 
particle  temperature,  T ^  the  bulk  temperature,  pp  the  apparent 
density  of  the  char,  Ag  the  specific  internal  surface  area  of  the  char 
particle,  and  /cr  is  the  intrinsic  reactivity.  The  effectiveness  factor  rj  is 
calculated  as 

*7  = -j2  (4>  coth  4>  -  1)  (4) 

4> 

where  </>  is  the  Thiele  modulus: 


<i P 


and  pox  is  the  density  of  oxidant  in  the  bulk  gas,  Sb  is  the  stoichiom¬ 
etric  coefficient  in  Reaction  (1)  in  kg02/kgchar  and  De  is  the  effective 
diffusion  coefficient  in  the  particle  pores.  Assuming  that  the  pore 
size  distribution  is  unimodal  and  the  bulk  and  Knudsen  diffusion 
proceed  in  parallel,  De  is  given  by 


6 


1  1 

+ 


n  -1 


D/cn  D 


0 


where  D0  is  the  bulk  molecular  diffusion  coefficient  and  6  is  the 
porosity  of  the  char  particle: 


9 


Pt 


where  Tp  is  the  particle  temperature,  rp  is  the  mean  pore  radius  of 
the  char  particle,  and  Mw  ox  is  the  oxygen  molecule  weight.  The  spe¬ 
cific  internal  surface  area  (Ag)  of  the  char  particle  is  assumed  in  this 
model  to  remain  constant  during  char  combustion.  The  intrinsic 
reactivity  (/q)  is  of  Arrhenius  form: 

hi  =  A;e-(£,'/i?Tp)  (9) 


where  A,  is  the  pre-exponential  factor  and  £,  is  the  activation 
energy. 

Now  we  are  facing  the  question  of  adapting  the  above  char 
combustion  model  to  the  Guasare  coal.  To  this  end  we  have  again 
used  the  IFRF  measurements  [23]  of  Guasare  char  burnout  as  a 
function  of  time  for  temperatures  of  950,  1300  and  1400  °C,  see 
Fig.  4(  right).  These  measurements  were  carried  out  at  4%  oxygen 
volume  fraction.  The  measured  char  morphology  data  were:  the 
apparent  char  density  equal  to  339  kg/m3  and  porosity  equal  to 
74%.  Fixing  the  mass  diffusion-limited  rate  constant  (Ci)  at  a  value 
of  5  x  10“12  rn3/K0  75  s,  the  char  porosity  at  the  74%,  and  the  spe¬ 
cific  surface  area  at  2.5  x  104  m2/kg  we  optimize  values  of  both  ki¬ 
netic  parameters  At  and  E{  to  obtain  a  proper  fit  to  the  IFRF  data,  as 
shown  in  Fig.  4(right).  The  optimization  results  in  the  values  given 
in  Table  5.  Thus,  in  our  CFD  predictions  presented  in  Paragraph  4, 
the  initial  particle  diameter  remains  unaltered  through  the  char 
combustion  processes  while  the  rate  of  char  oxidation  is  calculated 
using  Eq.  (3).  The  char  density  is  changing  during  the  char  burnout 
process  according  to  the  following  equation: 


B 


(1  -  do) 


where  p0  is  the  initial  density  of  char  particle,  a0  is  the  initial  ash 
content  and  B  is  the  burnout.  During  the  devolatilization  process, 
macromolecules  break  up  and  formation  of  tar  and  gases  can  lead 
to  swelling  of  the  particle.  In  our  model  the  swelling  coefficient 
(Csw)  takes  a  value  of  2. 

To  complete  the  model  description  we  quote  the  particle  energy 
balance  which  considers  both  the  convective  heat  transfer  and  the 
absorption/emission  of  radiation  at  the  particle  surface: 
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Table  5 

The  parameters  for  the  intrinsic  char  combustion  model  of  Guasare  coal. 


Parameter 

Symbol 

Value 

Unit 

Mass  diffusion-limited  rate  constant 

Ci 

5  x  1(T12 

m3  /K0-75  s 

Pre-exponential  factor 

A- 

1  x  1(T3 

kg/m2  s 

Activation  energy 

Ei 

5  x  107 

J/kmol 

Char  porosity 

0 

0.74 

- 

Mean  pore  radius 

Ep 

1  x  10~7 

m 

Specific  internal  surface  area 

Ag 

2.5  x  104 

m2/kg 

Tortuosity 

T 

C2 

- 

mpcp  =  hAp (T„  -  Tp )  +  epApa (g*  -  7^)  (11) 

where  mp  is  the  mass  of  the  particle  (kg),  cp  the  heat  capacity  of  the 
particle  (J/kg  K),  Ap  the  surface  area  of  the  particle  (m2),  Tp  the  par¬ 
ticle  temperature  (K),  T ^  the  local  temperature  of  the  continuous 
phase  (K),  h  the  convective  heat  transfer  coefficient  (W/m2  K),  sp 
the  particle  emissivity  (sp  =  0.8),  a  the  Stefan-Boltzmann  constant 
(W/m2  K4),  and  0R  is  the  radiation  temperature,  (K).  This  equation 
assumes  that  there  is  negligible  internal  resistance  to  heat  transfer, 
i.e.,  the  particle  is  at  a  uniform  temperature  throughout.  During  the 
char  combustion  the  additional  termfh  ■  ^  •  Hreac  has  to  be  added  to 
the  right  hand  side  of  Eq.  (1 1 ).  This  term  accounts  for  the  energy  re¬ 
lease  in  the  surface  reactions  ( Hreac )  which  is  directly  absorbed  by 
the  char  particle.  The  coefficient  fh  takes  a  value  of  0.3. 

In  summary  of  this  paragraph  we  refer  to  Fig.  5  which  shows  the 
coal  combustion  models  used  in  our  simulations.  The  devolatiliza¬ 
tion  is  modeled  using  the  CPD  model  with  the  model  parameters 
listed  in  Table  4.  The  volatile  matter  combustion  is  calculated  using 
Eqs.  (12)-(15)  while  for  the  char  combustion  rate  we  use  Eq.  (3) 
with  the  parameters  listed  in  Table  5. 

3.2.  Turbulence  and  chemistry -turbulence  interaction 

The  k-s  model  [30]  with  the  following  model  constants: 
Cie  =  1.44,  C2s  =  1.92,  Cu  =  0.09,  ok  =  1.0,  (7S  =  1.3  has  been  used 
for  the  flow  predictions. 

Combustion  of  the  volatiles  has  been  simplistically  represented 
by  two  overall  reactions: 

Ci.2H4.48O0.44  +  1.5O2  =  1.2CO  +  2.24H20  (12) 

CO  +  0.5O2  =  C02  (13) 

whose  rates  are  calculated  using  the  Eddy  Break-Up  (EBU)  model 
[31].  The  rate  of  production  or  consumption  of  species  i  is  then  gi¬ 
ven  by  the  smallest  of  the  two  expressions  below: 


Fig.  5.  The  coal  combustion  model  used  in  this  work. 


where  A  and  B  are  two  model  empirical  constants  equal  to  4  and  0.5, 
respectively.  YP  is  the  mass  fraction  of  any  product  species  (P)  and 
Yr  is  the  mass  fraction  of  the  particular  reactant  ( R );  v'i  r  is  the  stoi¬ 
chiometric  coefficient  for  species  i  in  reaction  r;  k  and  s  are  the  tur¬ 
bulent  kinetic  energy  and  the  turbulent  dissipation  derived  from 
the  k-s  model. 

3.3.  Particles  tracking  and  radiation 

In  discrete  phase  modeling,  coal  particles  of  known  size  distri¬ 
butions  and  properties  are  injected  into  the  combustion  chamber 
and  tracked  in  a  Lagrangian  fashion  throughout  the  computational 
domain.  The  trajectory  of  a  discrete  phase  particle  is  computed  by 
integrating  the  forces  acting  on  the  particle.  Only  the  drag  and 
gravity  forces  are  considered.  The  particle  tracking  procedure  used 
in  this  work  was  developed  at  the  beginning  of  the  eighties  at  the 
University  of  Sheffield  [32,33]  and  later  on  it  was  implemented 
into  the  FLUENT  code. 

The  Discrete  Ordinates  (DO)  [34,35]  radiation  model  is  used  to 
solve  the  radiative  heat  transfer  equation.  The  DO  radiation  model 
solves  the  Radiative  Transfer  Equation  (RTE)  for  a  finite  number  of 
discrete  solid  angles  each  associated  with  a  vector  direction 
fixed  in  the  global  Cartesian  system  (x,y,z).  The  DO  model  solves 
for  as  many  transport  equations  as  there  are  directions  Y.  The 
solution  method  is  identical  to  that  used  for  the  fluid  flow  and  en¬ 
ergy  equations.  The  DO  model  considers  the  radiative  transfer 
equation  in  the  direction  ~s  as  a  field  equation.  This  equation  is 
written  as 

V-  (/(  r  ,  s  )  s  J  +a/(  r  ,  s  )  =  a— — 1-^  /  /(  r ,  s')<P(  s  •  s')dQ' 

(16) 

where  Y  is  a  position  vector,  ~s  is  a  direction  vector,  s'  is  a  scatter¬ 
ing  direction  vector;  a  is  the  Stefan-Boltzmann  constant,  T  is  the  lo¬ 
cal  temperature,  I  is  radiation  intensity,  (P-phase  function  and  Q'  is 
the  solid  angle.  The  absorption  coefficient  (a)  is  prescribed  a  value 
of  1.5  m-1  and  kept  constant  through  the  furnace  volume.  The  value 
has  been  derived  using  the  narrow  angle  radiance  measured  by  the 
IFRF  [6],  see  Section  4.  The  scattering  has  been  omitted.  When  the 
solution  for  the  RTE  is  known  the  energy  source  term  for  energy 
equation  is  calculated  using  the  radiation  flux  qr: 

Srad  =  -Vqr  (17) 

3.4.  Nitric  oxides 

In  this  work  the  fuel,  thermal,  prompt  and  N20  paths  as  well  as 
NO  reburning  have  been  considered,  as  shown  in  Fig.  6.  The  trans¬ 
port  equations  for  nitric  oxide  (NO)  and  for  intermediate  species 
(HCN)  have  been  solved  for  while  the  N20  and  N-radical  concentra¬ 
tions  are  calculated  using  the  partial  equilibrium  assumption  (see 
below).  The  sources  for  HCN  and  NO  have  been  calculated  averag¬ 
ing  the  instantaneous  source  term  (S,  =f(T))  over  the  temperature 
fluctuations: 

rFmax 

(Si)  =  /  P(T)Si(T)dT  (18) 

using  a  PDF-/?  function.  The  single-variable  probability  density  func¬ 
tion  (PDF)  approach  used  in  our  work  was  originally  proposed  by 
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Hand  et  al.  [36],  extended  by  Peters  and  Weber  [37]  and  finally 
implemented  into  the  FLUENT  code.  T0  is  the  reference  temperature 
of  the  system,  Tmax  is  the  adiabatic  temperature  of  the  flame. 

Since  the  amount  of  NO  formed  is  small,  the  reactions  involved 
in  the  NO  chemistry  have  been  separated  from  the  pulverized  coal 
oxidation  mechanism.  Thus,  chemical  formation  and  reduction 
rates  of  NO  and  NO-concentrations  are  post-processed  using  previ¬ 
ously  computed  velocity,  turbulence,  temperature  and  chemistry 
fields.  As  shown  in  Fig.  6,  the  post-processing  involves  solving 
the  NO  and  the  HCN  balance  (transport)  equations  whose  source 
terms  are  calculated  as  described  below. 


3.4.2.  Fuel  path 

It  is  assumed  in  this  work  that  all  fuel  nitrogen,  both  from  vol¬ 
atiles  and  char,  is  converted  into  HCN  that  can  undergo  reduction 
or  formation  of  NO,  as  shown  in  Fig.  6. 

coal-N  -►  HCN  (19) 

HCN  +  02  -►  NO  H -  (20) 

HCN  +  NO  -►  N2  +  ---  (21) 


The  HCN  release  rates  are  given  by  Eqs.  (22)  and  (23),  where  mchar 
and  mVoi  are  the  combustion  rates  of  char  and  volatile  matter  (in 
kg/s),  respectively;  mchar  is  delivered  by  the  char  combustion  model 
(the  intrinsic  model)  and  mvo[  is  calculated  using  the  CPD  devolatil¬ 
ization  model;  V  is  the  cell  volume. 


r  ftlchar  '  WN ,char  '  -^HCN  .  /  3 

•JHCN,  char  = - iT7 - 77 -  (Kg/m  S) 


>HCN  ,vol 


MN-V 

TTlpol  '  tVjsj  ,vol  Mhcn 

mMv 


(kg/m3  s) 


(22) 

(23) 


The  HCN  depletion  rates  are  given  by  De  Soete  [38]  expressions: 

IT  =  k4XHcNXo2  exp  f-fif)  -jfi  (mol/m3  s)  (24) 

r2  =  MucnXn,  exp  ^  (mol/m3  s)  (25) 

where  k4  =  3.5  x  1010  s_1,  k5  =  3.0s_1,  E4  =  280.5  kj/mol, 
E5  =  251.2  kj/mol.  The  instantaneous  source  term  SHcn  reads  then 

^HCN  =  ShcN ’v°l  "*■  ^HCN, char  ~  (rl  +  r2)  ‘  -^HCN  10  3  (l<g/m3  S)  (26) 


and  is  averaged  over  the  temperature  fluctuations  using  Eq.  (18). 


The  NO  source  term  from  fuel  nitrogen  is  expressed  as  follows: 

Sm.fuei  =  (IT  -  r2)  •  Mn0  •  1 CT3  (kg/m3  s)  (27) 


3.4.2.  Prompt  path 

De  Soete  [38]  proposed  a  roughly  estimated  chemical  reaction 
rate  for  prompt  NO  formation: 

i6  =  /  -71  •  [02]6 [N2] [CxHy]  exp  '  V  (mol/m3  s)  (28) 

where  the  exponent  b  may  vary  between  0  and  1,  and  depends  on 
the  conditions  in  the  flame,  namely  the  local  mole  fraction  of  oxy¬ 
gen  [02]: 

'1.0,  [02]  <4.1  X  10~3 

I  —3.95  —  0.9  ln[02],  4.1  x  10~3  C  [02]  <  1.11  x  10^2 
]  -0.35  -  0.1  ln[02],  1.11  x  10~2  <  [02]  <  0.03 
„  0,  [02]  ^  0.03 

The  constant  A  and  activation  energy  Ea  take  the  following  values: 
A  =  6.4  x  \0\f)b+\  Ea  =  303.5  kj/mol.  The  prompt  factor  /  ac¬ 
counts  for  the  type  of  fuel  and  is  calculated  in  the  following  way 
[39]: 

/  =  4.75  +  0.0819  c- 23.2| +  32 T_  12.2 T 

2  X2  A3 

where  c  is  the  number  of  carbon  atoms  per  molecule  of  fuel  and  is 
assumed  to  be  c  =  1.2  for  the  volatile  matter  of  the  Guasare  coal 
and  X  is  the  air  excess  ratio. 

The  NO  source  term  for  prompt  path  reads  then 

•SnO,  prompt  =  ivMNo-10-3  (kg/m3  s)  (29) 

3.4.3.  Thermal  path 

The  extended  Zeldovich  mechanism  incorporates  the  following 


reactions: 

o  +  n2  k ^ 

NO  +  N 

(30) 

N  +  o2 

NO  +  O 

(31) 

N  +  OH 

NO  +  H 

(32) 

!  Does  not  appear  explicitly 
|  in  the  post-processor 


Fig.  6.  Path  of  NO  formation  and  reburning. 
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where  k+\ ,  k+2 ,  k+3  stand  for  the  forward  and  /<_! ,  /<_2 ,  /<-3  for  the 
backward  rate  constants  corresponding  to  Reactions  (30)-(32).  Un¬ 
der  the  assumption  of  partial  equilibrium  for  N-radical,  the  NO  for¬ 
mation  rate  can  then  be  calculated  as 

1  _  /<„i/<-2[NO]2 

r7  =  2 fc+,  [0]  [N2]  - - •  V  (mol/m3  s)  (33) 

1  +  /<+2  [02]+fc+3  [OH] 

The  0  and  OH  radical  concentrations  are  calculated  by  assuming 
equilibrium  of  the  following  reactions: 

02  +  H  ^  OH  +  O  (34) 

02  +  M  ^  M  +  0  +  0  (35) 

Finally,  the  NO  source  term  for  thermal  path  reads 

Sno, thermal  =  r7  ■  MN0  •  10“3  (kg/m3  s)  (36) 


3  A  A.  N20  path 

In  low  temperature  combustion  (below  1500  I<)  of  lean  mix¬ 
tures  of  hydrocarbons  (or  at  pressures  higher  than  the  atmo¬ 
spheric)  a  mechanism,  in  which  the  N20  molecule  is  the 
intermediary,  becomes  important  [40].  The  concentration  of  N20 
is  assumed  to  be  in  a  quasi-steady  state: 

IN  nl  =  Mn2][Q][M]  +  mno]2 

1  2  1  (kn2  +  k„3)[0]  +  (kn4  +  fc„5)[H]  1  ; 

Assuming  that  the  N20  is  converted  to  NO  the  rate  of  NO  formation 
is 

r8  =  2 [N20]  (/<„2  [0]  +  kn4 [H] )  •  V  (mol/m3  s)  (38) 

where 

[M]  =  1.4-  [02]  +  3.0  •  [C02]  +  1.7-  [N2]  +  12-  [H20] 


+  measurements 
—  simulation 


+  measurements 
—  simulation 


Fig.  7.  Velocity  (left)  and  temperature  (right)  profiles  along  the  measurement  traverses.  The  measured  values  taken  from  [6].  For  location  of  the  measurement  traverses  see 
Fig.  1. 


N.  Schaffel  et  al. / Combustion  and  Flame  156  (2009)  1771-1784 


1779 


Then,  the  NO  source  term  for  N20  path  is  calculated  as 

Sno,n2o  =  r8  •  Mno  •  10  3  (kg/m3  s)  (39) 

3.4.5.  Rebuming  in  gaseous  phase 

Following  the  work  of  Chen  et  al.  [41,42]  the  major  reburning 


reactions  are 

NO  +  CH2  ^  HCN  +  OH  (40) 

NO  +  CH  ^  HCN  +  O  (41) 

NO  +  C  ^  CN  +  O  (42) 

and  in  this  work  it  is  assumed  that  the  volatile  matter  is  the  reburn¬ 
ing  agent,  so  the  global  NO  reduction  rates  are 

r4  =  4  x  1CT4  •  (k9  +  ^10)  [Ci.2H4  4O0.44]  [NO]  (mol/m3  s)  (43) 
r5  =  4  x  1CT4  •  ku  [C,  2  H4  4O0.44]  [NO]  (mol/m3  s)  (44) 


where  fc9  =  5.30xl09-T“154-exp(-^§P),  k10  =  3.31  xl013T'3  33- 
exp (-1^9°),  k|  1  =  1.356  x  1012  •r~2,64-exp(-144'°90).  Then,  the  source 
term  for  NO  reburning  is  calculated  as 

SNO,(,om.re/,un,  =  -(r4  +  r5)-MNo-10^  (kg/m3  s)  (45) 

3.5.  Reburning  on  char  particles 

The  heterogeneous  reduction  of  NO  on  the  char  surface  is  mod¬ 
eled  following  the  work  of  Levy  et  al.  [43]: 

~itf)'Cs'Ag  (mol/m3  s)  (46) 

where  Ag  is  the  (BET)  specific  surface  area  and  cs  is  the  solid  matter 
concentration  and:  k6  =  2.27  x  10“3  mol/Pa  m2,  E6  =  142.7  kj/mol. 
The  source  term  for  NO-reburning  on  the  char  particles  is  expressed 
as  follows: 

SN0.hetjebwn  =  r3  •  MN0  •  10“3  (kg/m3  S)  (47) 

Thus,  considering  the  above-mentioned  mechanisms  of  NO  for¬ 
mation  and  destruction,  the  instantaneous  source  term  for  the  NO 
transport  equation  is  evaluated  as 

■SnO  —  ^NO  fuel  H-  ^NO, prompt  "T  ^NO, thermal  H-  -SnO,N20  +  ^NO, horn. reburn 

T-  ^NO het.reburn  (48) 

Obviously,  the  above  source  term  is  then  averaged  over  tempera¬ 
ture  fluctuations  using  Eq.  (18). 

4.  Results  and  discussion 

A  quarter  of  the  furnace  has  been  modeled  only  taking  advan¬ 
tage  of  the  symmetry.  An  unstructured  mesh  consisting  of 
700,000  cells  has  been  used  with  the  small  cell  size  being  around 
2  mm  and  the  largest  around  75  mm.  For  numerical  representation 
of  the  central  jet  inlet  area  600  cells  have  been  used  while  the 
small  fuel  jet  inlet  area  has  been  represented  by  150  cells.  Conver- 
gency  has  been  achieved  after  performing  around  50,000  itera¬ 
tions.  A  grid  sensitivity  study  was  performed  and  two  grids 
consisting  of  700,000  and  2,500,000  cells  were  tested.  Predictions 
using  both  grids  did  not  differ  by  more  than  few  percent  and  the 
smaller  grid  has  been  used  for  all  simulations  presented  below. 

4.1.  Flow  field  and  recirculation 

In  Fig.  7(left)  the  measured  and  calculated  velocities  along  the 
seven  traverses  are  shown.  At  the  first  traverse  the  computed  val¬ 


ues  correspond  well  with  the  measured  velocities  except  for  the 
comburent  centerline  velocities  which  were  measured  about 
20  m/s  lower  than  the  calculated  values.  Here  the  measured  values 
are  obviously  incorrect.  For  a  128  mm  diameter  jet  of  65  m/s  inlet 
velocity,  at  x/d  =  1.2  which  is  in  the  jet  potential  core,  the  velocity 
must  be  still  around  65  m/s.  At  the  first  two  traverses  the  combu¬ 
rent  jet  can  be  clearly  distinct  from  the  coal  jet.  The  coal  jet  merges 
into  the  comburent  jet  near  the  third  traverse.  Downstream  of  the 
third  traverse  both  jets  merge  and  form  one  stream  that  can  be  ob¬ 
served  at  the  next  traverses.  At  a  distance  of  2  m  downstream  of 
the  furnace  front  wall  (from  the  fifth  traverse  onwards)  the  veloc¬ 
ity  profile  is  flat.  As  shown  in  Fig.  7(left)  the  size  of  the  recircula¬ 
tion  region  formed  in  the  furnace  and  its  low  magnitude 
(negative)  velocities  are  well  reproduced  in  the  computation.  Gen- 


r3  =  -k6  •  XN0  •  p  ■  exp 
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erally,  the  velocity  predictions  (the  flow  field  together  with  the 
recirculation)  are  in  good  agreement  with  the  measurements. 

4.2.  Temperature  field  and  radiative  heat  fluxes 

The  temperature  profiles  along  the  traverses  are  given  in 
Fig.  7(right).  No  substantial  difference  between  the  model  predic¬ 
tions  and  the  measurements  can  be  seen  with  exception  of  the  fuel 
jet  at  the  first  traverse.  The  temperature  level,  as  well  as  the  temper¬ 
ature  peaks  are  reproduced  well  by  the  simulations.  Ignition  of  the 
fuel  jet  takes  place  somewhere  between  the  first  and  the  second  tra¬ 
verse  so  that  the  peak  temperature  of  around  1 800  I<  occurs  between 
the  second  and  third  traverses.  This  feature  is  again  well  reproduced 
in  the  numerical  simulations.  From  the  fourth  traverse  downwards, 
the  temperature  profile  is  flat  since  slow  combustion  proceeds 
downstream  in  the  furnace  and  the  temperature  level  is  everywhere 
around  1600  K.  Generally,  the  agreement  between  the  temperature 
predictions  and  the  measurements  is  quite  satisfactory. 


To  obtain  such  a  good  accuracy  of  the  temperature  predictions 
it  is  imperative  that  the  value  of  the  absorption  coefficient  appear¬ 
ing  in  Eq.  (16)  is  properly  selected.  Following  the  method  of  Lalle- 
mant  et  al.  [44]  and  Sayre  et  al.  [45]  we  use  here  the  narrow  angle 
radiance  measured  at  the  fifth  traverse  (2.05  m).  Fig.  8  shows  that 
using  a  value  of  1.5  m-1  which  accounts  for  the  C02/H20  radiation 
as  well  as  for  soot  and  particles  radiation,  the  predicted  and  mea¬ 
sured  radiation  intensities  are  in  good  agreement.  The  verification 
of  the  correctness  of  the  radiative  heat  transfer  procedure  is  dem¬ 
onstrated  in  Fig.  9  showing  very  good  agreement  between  the 
measured  and  the  predicted  total  incident  heat  flux  at  the  furnace 
wall.  Flatness  and  high  values  of  the  radiative  heat  flux  are  charac¬ 
teristics  of  mild  combustion  [46,47]. 

4.3.  Oxygen  and  carbon  dioxide  concentration  fields 

In  Fig.  10(left)  the  profiles  of  oxygen  concentration  are  shown. 
No  major  differences  between  the  measurements  and  calculated 
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Fig.  10.  Oxygen  (left)  and  carbon  dioxide  (right)  concentration  profiles  along  the  measurement  traverses.  The  measured  values  taken  from  [6].  For  location  of  the 
measurement  traverses  see  Fig.  1. 
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values  can  be  observed.  However  the  oxygen  concentrations  near 
the  coal  nozzle  are  too  low  in  comparison  with  the  experimental 
values  which  is  consistent  with  the  temperature  predictions  in  this 
region,  see  Fig.  7(right).  It  is  interesting  to  note  that  the  same  mod¬ 
el  deficiency  was  already  observed  in  numerical  simulations  of  the 
natural  gas  combustion  [47]  and  it  was  related  to  inaccurate  pre¬ 
dictions  of  the  entrainment  of  the  fuel  jet  [48].  Downstream  of 
the  fourth  traverse  the  oxygen  concentration  is  uniform  and  prac¬ 
tically  everywhere  in  the  furnace  is  about  3%.  Altogether,  the 
numerical  predictions  of  oxygen  concentration  and  the  measure¬ 
ments  are  in  good  agreement.  Carbon  dioxide  profiles  are  given 
in  Fig.  lO(right).  The  calculated  values  are  very  close  to  the  mea¬ 
sured  values  along  all  traverses. 

4.4.  Carbon  monoxide  concentration  field 

Carbon  monoxide  profiles  are  illustrated  in  Fig.  11  (left).  In  our 
simulations  carbon  monoxide  is  generated  from  the  combustion 
of  volatiles  only.  If  oxygen  is  available  it  can  be  further  oxidized 
to  carbon  dioxide.  It  is  formed  mostly  within  the  coal  jet  up  to 
the  forth  traverse.  The  highest  concentration  of  CO  is  about  6% 
and  it  is  located  at  the  third  traverse.  The  prediction  of  its  peak  va¬ 
lue  agrees  perfectly  with  the  measurement.  Downstream  of  the 
fifth  traverse  no  considerable  amount  of  carbon  monoxide  is 
detectable.  Although  the  carbon  monoxide  predictions  are  in  very 
good  agreement  with  the  measured  values,  we  regard  this  agree¬ 
ment  as  coincidental  for  several  reasons.  The  monoxide  is  formed 
and  oxidized  in  a  number  of  sequential  and  parallel  reactions.  It 
is  released  as  a  component  of  the  volatiles,  as  a  product  of  the  vol¬ 
atiles  oxidation,  and  as  a  product  of  char  combustion.  The  monox¬ 
ide  is  oxidized  by  OH  radicals  as  well  as  by  O  radicals,  although 
with  much  lower  a  rate.  Furthermore,  the  water  gas  shift  reaction 
alters  CO  concentrations.  None  of  these  processes  are  taken  into 
account  in  our  modeling.  Thus,  from  the  fundamental  point  of 
view,  the  works  of  Kim  et  al.  [18,19]  are  superior  since  their  at¬ 
tempt  to  take  these  CO-involving  reactions  into  account  although 
the  CO-predictions  in  Kim’s  simulations  depart  from  the  measured 
values. 

4.5.  Nitric  oxide  concentration  field 


see  Fig.  12(left),  one  can  observe  that  the  NO  reburning  rate  seems 
to  be  over-estimated.  Here  we  recall  that  in  our  model  the  entire 
volatile  matter  is  involved  in  the  NO-reburning  (see  Eqs.  (43)  and 
(44))  which  may  be  questionable.  The  second  reason  is  related  to 
the  lack  of  oxygen  in  the  fuel  jet  obtained  in  our  predictions  while 
in  the  measurements  still  around  2%  of  oxygen  can  be  observed. 
The  oxygen  concentration  affects  strongly  the  NO-fuel  formation 
path  (see  Eq.  (24)).  For  the  sake  of  completeness  we  report  that 
the  thermal-NO  source  (Eq.  (36)),  the  prompt-NO  source  (Eq. 
(29))  and  the  N20-N0  source  (Eq.  (39))  are  in  the  range  of 
0 6.6  x  10“8,  0-f1.6x1(T6,  0  -f  4.6  x  10-9  kmol  m3  s,  respec¬ 
tively.  Integration  of  the  NO  sources  over  an  entire  furnace  volume 
provides  a  further  insight  into  the  NO-mechanisms.  The  overall  NO 
production  rate  through  the  fuel  mechanism  is  calculated  to  be 

18.5  x  10  6  kgN0/s  (see  Table  6).  The  NO  reburning  rate  is 

11.5  x  10  6  kgN0/s  so  that  the  net  NO  production  rate  in  the  fur¬ 
nace  amounts  to  7  x  10  6  kgN0/s.  The  difference  between  the 
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Nitric  oxide  concentration  profiles  are  given  in  Fig.  12(left) 
showing  that  the  concentrations  of  NO  are  reasonably  well  pre¬ 
dicted  by  the  model,  not  only  qualitatively  but  also  quantitatively. 
It  can  be  observed  that  the  NO  formation  begins  upstream  of  the 
first  traverse,  however  most  of  the  NO  is  formed  upstream  of  the 
third  traverse  in  the  volatiles  released  zone.  At  the  second  traverse 
the  highest  peak  of  900  ppm  has  been  measured.  This  peak  is  also 
well  reproduced  in  the  calculations.  Downstream  of  the  fourth  tra¬ 
verse  the  nitric  oxide  profile  is  flat  and  of  a  low  level  (around 
300  ppm).  At  the  outlet  333  ppm  NO  was  calculated  and  this  value 
agrees  very  well  with  the  measurements  (320  ppm). 

Among  the  six  terms  which  appear  in  Eq.  (48)  for  the  net  NO 
formation  rate  only  the  fuel-NO  source  (Eq.  (27)),  the  homogenous 
NO-reburning  source  (Eq.  (45))  and  heterogenous  NO-reburning 
source  (Eq.  (47))  are  significant  and  there  are  plotted  in  Fig.  12. 
The  fuel-NO  source,  SNOfUei,  is  of  the  order  of  0  -f  7.8  x  10“4  kmol/ 
m3  s  while  the  homogenous  NO-reburning  rate,  Sm_hom  reburni  is  in 
the  range  of  -2.7  x  10“4  -f  0  kmol/m3  s  and  the  heterogenous 
NO-reburning  rate,  Sm_hetreburn,  is  in  the  range  of  -2.0  x  1CT4-f 
0  kmol/m3  s.  Thus,  these  three  sources  are  of  the  same  order  of 
magnitude  and  the  in-flame  NO-concentrations  are  determined 
by  their  balance  which  is  then  altered  due  to  the  NO  advection 
and  diffusion.  The  NO  post-processor  predicts  a  rapid  reburning 
occurring  already  in  the  fuel  jet  where  NO  is  formed  directly  from 
HCN  decomposition.  Examining  the  NO  predictions  in  the  fuel  jet, 


Fig.  11.  Carbon  monoxide  concentration  profiles  along  the  measurement  traverses. 
For  location  of  the  measurement  traverses  see  Fig.  1. 
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Fig.  12.  Nitric  oxide  (left)  along  the  measurements  traverses.  Sources  in  the  NO  balance  equation  (right).  From  the  top  to  the  bottom:  SN0 jue/  -  Eq.  (27),  SN0,het.rebum  -  Eq.  (47), 

S NO.hom. reburn  ~  Eq.  (45). 


Table  6 

Mass  balance  of  NO. 

MN0,in  (kgwo/s)  4.0  X  10'6 

Mm.out  (kgN0/s)  -1 1 .0  x  10'6 

M NO. fuel  (kgN0/s)  18.5  X  10-6 

Mmsebum  (kgN0/s)  -11.5  X  10“6 


amount  of  the  NO  at  the  furnace  exit  (11  x  1CT6  kgN0/s)  and  at  the 
furnace  inlet  (4  x  1CT6  kgN0/s)  amounts  exactly  to  7  x  1CT6  kgN0/s. 
Thus,  only  two  terms  in  the  overall  NO-balance  are  significant;  the 
NO  generation  via  the  fuel  mechanism  and  the  NO  reburning 
mechanism. 

4.6.  Char  burnout 

Fig.  13  shows  both  the  measured  and  predicted  char  burnout 
along  the  centerline  of  the  fuel  jet.  The  char  burnout  predictions 
are  in  good  agreement  with  the  measured  data  up  to  80%  burnout. 
For  higher  degrees  of  burnout  the  model  over-predicts  the  char 
oxidation  rates  and  at  the  furnace  outlet  a  complete  burnout  is  pre¬ 
dicted  (see  Table  7).  This  is  a  consequence  of  the  fact  that  the 
parameters  of  the  char  combustion  model  have  been  derived  using 
the  measured  data  (Fig.  4( right))  that  extend  up  to  85%  burnout.  As 
can  be  seen  in  Fig.  4(right),  an  extrapolation  of  the  char  model  pre¬ 


dictions  to  residence  times  of  5-6  s  results  in  complete  burnout.  In 
order  to  predict  the  last  stages  of  burnout,  corrections  to  the  char 
model  are  needed  to  slow  down  the  rate  as  the  char  oxidation  pro¬ 
ceeds  (see  for  example  [49,50]). 

4.7.  Furnace  outlet 

The  temperature,  oxygen,  carbon  oxides  and  nitric  oxide  con¬ 
centrations  and  char  burnout  at  the  outlet  of  the  furnace  for  calcu¬ 
lations  and  measurements  are  summarized  in  Table  7.  The  furnace 
exit  temperature  is  very  close  to  the  measured  value  and  the  differ¬ 
ence  is  about  50  K  which  corresponds  to  2.37%  of  the  furnace  ther¬ 
mal  input. 

5.  Conclusions 

CFD-based  numerical  simulations  of  the  IFRF  experiment  on 
MILD  combustion  of  coal  have  been  carried  out.  The  burner  system 
consists  of  a  central  jet  of  vitiated  air  preheated  to  1623  K  and  two 
jets  of  pulverized  coal.  This  is  indeed  a  unique  burner  design  as  for 
coal  combustion.  Prior  to  performing  the  numerical  simulations 
substantial  efforts  have  been  allocated  to  an  accurate  modelling 
of  combustion  of  Guasare  coal  which  has  been  used  in  the  IFRF 
experiment.  To  this  end  both  the  Chemical  Percolation  Devolatil- 


N.  Schaffel  et  al. / Combustion  and  Flame  156  (2009)  1771-1784 


1783 


Fig.  13.  Char  burnout  (left)  and  carbon  in  ash  (right)  along  the  centerline  of  the  fuel  jet.  The  measured  values  taken  from  [6]. 


Table  7 

Computed  and  measured  values  at  the  furnace  exit. 


Parameter 

Unit 

Experiment  [6]  Predictions  (this  work) 

T 

K 

1503 

1555 

o2 

vol.%  dry 

3.1 

2.52 

CO 

vol.ppm  dry 

<50 

10 

co2 

vol.%  dry 

25.50 

23.93 

NO* 

vol.ppm  dry 

320 

333 

Char  burnout 

o/ 

Zo 

99.4 

100 

Carbon  in  ash 

O/ 

Zo 

14.95 

None 

ization  (CPD)  model  and  the  char  combustion  intrinsic  reactivity 
model  have  been  adapted  to  Guasare  coal  combusted.  To  ensure 
accurate  predictions  of  radiative  heat  transfer,  an  appropriate  va¬ 
lue  of  the  local  absorption  coefficient  appearing  in  the  Radiative 
Transfer  Equation  has  been  derived  using  the  measured  data.  The 
detailed  validation  of  the  predictions  against  the  IFRF  in-furnace 
measurements  has  led  to  the  following  conclusions: 

•  The  flow-field  as  well  as  the  temperature  and  oxygen  fields  have 
been  accurately  predicted  by  the  CFD-based  model.  The  pre¬ 
dicted  temperature  and  gas  composition  fields  have  been  uni¬ 
form  demonstrating  that  slow  combustion  occurs  in  the  entire 
furnace  volume.  Thus,  this  realization  of  MILD  combustion  is 
inherently  different  to  fast  combustion  in  high  gradient  fields 
generated  in  swirl-stabilized  flames. 

•  The  CFD-based  predictions  have  highlighted  the  NOx  reduction 
potential  of  MILD  combustion  for  N-containing  fuels  through 
the  following  mechanism.  Before  the  coal  devolatilization  pro¬ 
ceeds,  the  coal  jet  entrains  a  substantial  amount  of  flue  gas  so 
that  its  oxygen  content  is  typically  not  more  than  3-5%.  The  vol¬ 
atiles  are  given  off  in  a  highly  sub-stoichiometic  environment 
and  their  N-containing  species  are  preferentially  converted  to 
molecular  nitrogen  rather  than  to  NO.  Furthermore,  there  exist 
a  strong  NO-reburning  mechanism  within  the  fuel  jet  and  in 
the  air  jet  downstream  of  the  position  where  these  two  jets 
merge.  In  other  words,  less  NO  is  formed  from  combustion  of 
volatiles  and  stronger  NO-reburning  mechanisms  exist  in  the 
MILD  combustion  if  compared  to  conventional  coal  combustion 
technology. 

•  Although  the  carbon  monoxide  predictions  happen  to  be  in  very 
good  agreement  with  the  measured  values  we  regard  this  agree¬ 
ment  as  coincidental  (see  text)  and  more  elaborated  schemes  of 


volatile  matter  combustion  are  desired.  The  char  burnout  figures 
have  been  well  predicted  until  80%  of  burnout.  However,  correc¬ 
tions  to  the  char  model  are  needed  to  slow  down  the  rate  as  the 
char  oxidation  proceeds  further. 
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